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O 

Q ' Abstract 

We study bosonisation in the massive Thirring and sine-Gordon models at finite 
temperature T and nonzero fermion chemical potential fi. For that purpose we use 
both canonical operator and path integral approaches, paying particular attention 
to the issues of thermal normal ordering and renormalisation. At T > and 
fj, = 0, the massive Thirring model bosonises to the sine-Gordon model with the 
same T = identification between coupling constants. We prove that not only the 
partition functions of the two models coincide, as was recently shown, but also that 
qq \ thermal averages of zero-charge operators can be identified. In particular, analysis 

ON ' of the point-split regularised fermion current then leads to the thermal equivalence 

between sine- Gordon kinks and Thirring fermions. At \x ^ 0, T > and working in 
Oh, perturbation theory about the massless Thirring model, we show that the bosonised 

Q_i| theory is the sine-Gordon model plus an additional topological term which accounts 

for the existence of net fermion charge excitations (the fermions or the kinks) in the 
thermal bath. This result generalises one recently obtained for the massless case, 
and it is the two-dimensional version of the low-energy QCD chiral Lagrangian at 
^ ■ finite baryon density. 



1 Introduction 

The sine- Gordon (SG) and massive Thirring (MT) models are just two of many different 
models which have been widely studied in two space-time dimensions (2D). They provide 
one of the simplest examples of bosonisation, in which a fermion theory (MT) is equiv- 
alent to a boson one (SG) with certain identifications between their coupling constants. 
Furthermore the weak limit of one theory is the strong limit of the other and conversely 
so that these models also give a simple example of duality, and the particle spectrum 
can be mapped from one model onto the other. This equivalence was first analysed in || 
using canonical operator methods and, later on, using path integral in [FJ, |3|]. In the MT- 
SG system, bosonisation takes place for the full spectrum. This should be contrasted, 

*E-mail: Gomez@eucmax.sim.ucm.es 
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for instance, to QCD where it happens only at low energies for which strongly coupled 
quarks and gluons are confined into hadrons. There the effective boson theory is the 
chiral Lagrangian for the lightest mesons which are the Nambu-Goldstone bosons of the 
chiral symmetry breaking pattern, and to lowest order it is the non-linear sigma model 

It is important to note that the particle spectrum of the bosonised theory may still 
contain fermion-like excitations. In the SG model, such excitations are soliton-like, cor- 
responding classically to nondissipative solutions of the equations of motion. One can 
assign a conserved charge to those configurations, corresponding to the fermion charge |5| 
which is a super-selection rule, so that the theory can be classified into different charge 
sectors. For instance, a possible one-fermion solution is the static kink |J. The particle 
spectrum of the SG theory is then constructed from the mesons (including kink-anti- 
kink 'breather solutions') and the kinks, which are in one to one correspondence with 
the Thirring model fermion-antifermion bound states and fermions, respectively. The 
situation is entirely analogous to low-energy QCD where the mesons are the pions, kaons 
and so on and the solitons are the skyrmions J7|, corresponding physically to the QCD 
baryons ||. An example in which the bosonised theory does not contain fermion modes 
is the Schwinger model (QED in 2D) which bosonises to a free massive scalar with mass 
proportional to the electric charge [[J. This scalar is a fermion-antifermion bound state 
in which the electric charge is confined. 

The link with QCD discussed above is just one of the motivations to study these 2D 
theories as toy models. However, one should recall that some theories in 2D also describe 
experimentally observed phenomena. Indeed, the SG equation is behind a device which 
produces some of the highest energy microwave signals available, having technological 
applications, for example, in space physics |TT], [TTJ . Such Josephson junctions consist of 



two layers of super-conducting materials separated by a very thin dielectric barrier, and 
Josephson tunnelling of the Cooper pairs across this barrier results in the relative phase 



angle of the two superconductors satisfying a SG equation Jl0|. As the two superconduc- 
tors are taken through the phase transition, kinks are formed in the junction and these 
are observed experimentally (T^|. The kinks are topological defects, the 2D analogues of 



cosmic strings or vortices, line-like defects formed in 4D when a system goes through a 
symmetry breaking phase transition (say of some group G to a subgroup H) for which the 
first homotopy group ni(G/H) ^ [14[ . In the early universe they may provide an 
explanation for the formation of structures such as galaxies [Q. An important problem 



in this context is to estimate the number density of such strings formed at the phase 
transition and also their length distribution. Recently progress has been made by using 
the analogy between experimentally observable systems such as 3 He and 4 He and the 
early universe fll. It would seem, however, that Josephson junctions provide an even 



more simple experimentally accessible system with which one could try to test ideas of 
defect formation. Furthermore the analysis should be simplified in this case by exploiting 
the duality between the SG and the MT models flEfl . For that purpose, one needs to 
know first what is the exact relation between these two models at T > and fi ^ 0, and 
how kink densities and fermion densities are related (see section |5] below). 

In this paper we therefore study bosonisation in the SG and MT models at finite 
temperature T and non-zero chemical potential /i. Recall that regarding bosonisation, 
temperature and chemical potentials play a completely different role. Temperature affects 
the boundary conditions of the propagator and hence bosonisation implies a change from 
Fermi-Dirac to Bose-Einstein statistics. However, a nonzero chemical potential allows 
the system to have net charged states. As we have commented before, this is especially 
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relevant when the bosonised theory admits fermion states as in the SG model or in chiral 
Lagrangians. In fact, we expect the chemical potential to reflect that property. Besides, 
H 7^ spoils the hermiticity of the Dirac operator and breaks charge conjugation (C) 
invariance. Indeed breaking the combination of C and Lorentz covariance (the thermal 
bath selects a preferred frame) can give rise to topological /i-dependent new terms in the 
bosonised action [17|, [18|, |19[| . 

We now introduce the models. In Minkowski space with metric (+, — ), the Lagrangian 
densities for the SG and massive Thirring models are respectively 

C MT $,?P}=i^-rno)4> + ^g 2 ]»(x)f(x), (1.2) 

where j M (x) = , ij;(x) r y tl i/;(x), a>o, m and 70 are bare parameters and and i[) are bosonic 
and fermionic fields respectively. The Minkowski gamma matrices are 

7o= ( 1 I )' 7l= ( -1 )' ^ = 7o7i= ( "q 1 J ) 



C SG [<f>} = + T§ cos M - 7o, (LI) 



so that 



{7^} = 2<r 7^75 = -e^7" 
with e i = 1. 

As we have commented before, the SG equations of motion admit nondissipative 
solutions for which the classical field <f) c i{t, ±00) = 2irn±/\, with n± £ 7L (with properly 
normalised 70). In fact, those classical static solutions with n + ^ are the kinks we 
discussed above. It is also worth noticing that the SG Lagrangian in (|1 . 1|) corresponds 
to a non-linear sigma model in 2D for a single Nambu-Goldstone-like field 0. Although 
there is no spontaneous symmetry breaking in 2D |2(J , the potential term in ( |1.1| ) breaks 



explicitly the symmetry — > + a with a £ IR (which we will call the chiral symmetry 
for reasons to become clear below) still preserving the symmetry — > + 27m/ A with 
n £ Z. These two symmetries are, respectively, the counterparts of the chiral and 
isospin symmetries for QCD, a and A playing the role of the pion mass squared and the 
inverse of the pion decay constant respectively. On the other hand, the chiral symmetry 
transformations in terms of the Thirring fermion are ip — > exp(ia\~f5)ip . The massless 
Thirring model is chirally invariant, the fermion mass term breaking that symmetry in 
the same way as the a term does in the SG Lagrangian. 

We now recall the results at T = fi = 0. In []]], Coleman, using a canonical field 
operator approach, considered chiral invariant combinations of zero-charge operators!] 
and showed that the SG and the massive Thirring models are equivalent in the following 
sense 

a 

mip(x)ip(x) = — cosA0(x), (1.3) 

3n{x) = —e^d u <p(x), (1.4) 

provided the constants of the two models are related through 

\ 2 1 

(1.5) 



47T 1 + g 2 / 71 "' 



a 
A 2 



pm, (1.6) 



1 These operators have the same number of fermions and anti-fermions in the MT, and correspond to 
operators which are invariant under (j> — > <fi + 2irn/\, with n e 2, in the SG model. 
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where p is the renormalisation scale and all the coupling constants and fields in ( |1.3|) -(|Ll) 



are renormalised (see sections |2.2| and |3.2| ). In fact the renormalisation convention chosen 
here is not exactly the same as that of Q (see below) and it also differs from by a 7r 
factor which may be absorbed in p — > p/ir in the above equation. All the physical results 
must be independent of the scale. Throughout this paper, we will refer to ( |1.3|) and ( |1.4[) 
as the mass and current equivalences respectively. It is very important to observe that 
whilst it may be inferred from Jl| that the above equivalences are strong, this is not the 
case — they are weak identities, only holding between the vacua of the two theories. We 
also recall that (|1.4j) allows one to establish that the net fermion charge (fermions minus 
anti-fermions) in the massive Thirring model corresponds to the net kink number (kinks 
minus anti-kinks) in the SG modelQ [[[], Equations (|1.5|) - (|L6| ) were also obtained using 
path integral techniques in |3|. In both |TJ and |3|, the approach followed was to 
expand formally around the «o — in the SG model and the mo = in the Thirring 
model — we will do the same here. An alternative approach would be, for instance, to 
work within a semiclassical approximation jlL |22| , which we will not attempt in this work. 



A question which may be asked is whether the equivalences ( p..3|) -(1.4) also hold in 



a heat bath where the vacuum expectation values are replaced by thermal expectation 
values. For bosons, 

< • > = i^rTr (•e~ p6 ) = f d<f> • exp(-5[0]), (1.7) 



Z(T) v ' Z{T) . ' periodic 

where Np is an infinite T-dependent constant arising in the path integral description fl23f , 
(3 = 1/T, the path integral boson fields are periodic in Euclidean time <f>(i + (3) — (f>{t) 
with i = it and 

Z{T) = Tr (e~^) = N fi f d(j> exp(-5[0]) 

^ ' J periodic 

is the partition function. Thermal averages and partition functions for fermions are 
defined analogously, but integrating over anti-periodic field configurations instead. It is 
important to remember that for p — 0, one works in the canonical ensemble, fixing the 
net fermion charge to zero. Consistently, in that case, the above trace is taken in the 
SG model over all possible states with zero net kink number or, equivalently, over field 
configurations that vanish at spatial infinity (configurations such as a kink-antikink pair 
are permitted). On the other hand, for p ^ 0, we are working in the grand-canonical 
ensemble and the net fermion number is not fixed. Hence the net averaged fermion density 
is nonzero and we are summing over all possible states in all the super- select ion sectors 
of the theory, i.e, over all possible net kink numbers in the SG model. In other words, 
the complete set of states over which the trace is taken can be chosen to be eigenstates 
of the Hamiltonian for p = and eigenstates of both the Hamiltonian and the number 
operator for p ^ 0. In |24|] path integral methods were used to show that the SG and 
MT model partition functions are equivalent at T > but p = 0. On the other hand, it 



was shown in [19| that for m = and p ^ 0, the Thirring model partition function is 
not the free boson one (oiq = 0), but it acquires an extra /i-dependent term accounting 
for the nonzero fermion net number in the thermal bath. 

In this work we analyse several aspects of the thermal bosonisation for this system 
which we believe have not been studied before. In the first part of the paper (section H) 
we extend the work of Coleman |l| to T > and p = 0, working with a canonical operator 



2 Note that here we are calling quantum configurations with any nonzero net charge and space-time 
dependence kinks, although strictly speaking the kinks are just time-independent classical configurations 
with net charge equal to one. 
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approach. Particular attention is paid to the question of renormalisation and the defini- 



tion of normal ordering in a heat bath. In the process, the results of |24| are reproduced 
though the approach is entirely different. The rest of the paper uses path integral meth- 
ods, which are introduced in section |3|. There we again prove that the partition functions 
of the two models are the same, but we pay special attention to renormalisation and to 
the connection with our results in the canonical formalism. The main objective of section 



is to prove that thermal averages of correlators of the zero-charge operators ( |1.3| )-( |lT4"D 



evaluated at different space-time points coincide. Such correlators cannot be inferred in 
general from the partition function which contains information about global thermody- 
namic observables such as the pressure or the condensates, but not about correlators, 
which physically yield for example thermal correlation lengths. In our case, the thermal 
version of ( |1.4| ) will constitute the thermal Thirring fermion-SG kink equivalence. This 
equivalence will appear again in section |5], where we analyse the case T > and \x 7^ 
for the partition function, thus extending [O] to the massive case and [El] to fi 7^ 0. The 



massless and fi = case was analysed in |25[ . A summary of the results presented in this 
paper may be found in [E6|. 



2 Bosonisation in the canonical operator approach 

at T > and fi = 

Here we work in Minkowski space with a canonical operator approach and extend the 



results of |T| to T > through the use of thermal normal ordering introduced in | 27| 
and defined below. Calculations are carried out in perturbation theory, expanding about 
cto = in the SG model and about mo = in the Thirring model. We first summarise 
useful results for 2D free partition functions and propagators at finite temperature and 
zero chemical potential. 

2.1 Results for free fields in 2D at T > and fi = 0. 
2.1.1 Notation and free partition functions 

For simplicity we take the spatial dimension to be finite < x < L with corresponding 
discrete momenta k := k n = 27m/ L. The infinite length limit is reached by substituting 
J2k — > L f-m f^- As opposed to the fields in the path integral (sections 0-0), field 
operators are entirely independent of temperature T = (3~ l and hence do not satisfy any 
periodicity conditions in imaginary timeQ. As usual, all the thermodynamic observables 
are obtained from the partition function and its derivatives by dividing by (3L. 

Consider a free real scalar field (f) of mass fi and a free massless fermion field ip with 
corresponding Hamiltonians 



M= ( L dx 
Jo 



7T 2 1 (d^ ^ 

2 2 [dx I 2 



f L 

; Hn = — i I dx 
Jo 



(2.1) 



where n = <9o0, V> = V^7o and the subscript denotes the free case. As usual, the free field 
operators (or more generally operators in the interaction picture) are written in terms of 

3 Temperature dependence, and thus questions of periodicity, only appear when considering thermal 
expectation values of operators. 
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annihilation and creation operators as 



tp(t,x) 



k 

1/2 

k 



-1/2 



ak,jxe 



-ik-x 



+ a 



ik-x 



(2.2) 



L- 1 / 2 \c k e- ik - x + d\e lk - x \ u{k) 



where uj\^ 



k + pr and the two-vector u{k) is 

9{-k) 
9(k) 



Also k ■ x is shorthand for kgt 




and 



a fc,A' a fc',/iJ — {cfc,4'} — {4,4'} — $k,k> 
with all other commutators and anti-commutators vanishing. Note that for the scalar 
field, our notation differs from the conventional one in that the annihilation and creation 
operators carry an extra label which is the mass p appearing in the free Hamiltonian. In 
terms of annihilation and creation operators, the boson Hamiltonian is as usual 



-"0 



£ 

k 



dk,p,d\.ji + a k pak,fi 



(2.3) 



where the standard normal ordering operation, : • which places the annihilation 
operators dk ^ to the right of creation operators a' k has been used to separate out the 
infinite vacuum energy. Note again the mass label on the normal ordering operation — it 
will be important below where objects such as : Hq : p with p ^ p are considered. (By this 



we mean that the operators must be expressed in terms of a mass p in ( |2.2| ), and then 
the creation and annihilation operators dk, p are ordered. This is equivalent to considering 
a mass shift in the Hamiltonian: H[f = Hq + (p 2 — p 2 )0 2 /2.) From ( |2.3| ), the free boson 
partition function is, in the limit L — + oo, 



log Z£ B {T) 



log [Tr{e~^o} 



—L 



-°° dk 
x> 2^ 



fan, 



+ log ( 1 - e 



(2.4) 



The first term in the far r.h.s above is divergent and corresponds to the vacuum energy. In 
the canonical approach it is renormalised by normal ordering. In path-integral methods 
there is no operator normal ordering and infinite vacuum terms are not removed, although 
they are irrelevant for the thermodynamics (they are T-independent) and can be ignored 
|23| . For the massless fermion hamiltonian one has 



whence 



£ 1*1 (44 - 44) = £ 1*1 (4 g fe + 44) - £ 1*1 

k 

+ log (l + e~^ k 



: H n 



£1*1, 



log 3 



2L / — 



00 dk 
2~ix 



(2.5) 



Note that since 



r+°° dk , / mhl \ r+°° dk , / 

/ —log(l-e-^l) =-2 / — log(l + e 



-P\k\ 



7TT 

IT' 



it follows that, ignoring vacuum terms, 



Z£(T) = Z° Q > B (T) = exp 



ttLT 
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(2-6) 



This result provides one of the simplest examples of bosonisation, and we will make use 
of it later on since it forms the basis of the equivalence of the SG and MT partition 
functions at T > 0. 
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2.1.2 Thermal normal ordering and thermal propagators 



Central to our discussion below will be the use of 'thermal normal ordering' 27]. The 
reason one is lead to consider a different normal ordering at T > comes from the fact 
that whilst at zero temperature, dealing with normal ordered products is easy as by 
construction their vacuum expectation value vanishes, for example, 



(0,p,\ : d Kfl d\^ : A |0,/2) 
this is no longer the case at T > where 



0. 



(2.7) 



N k>ii + 0. (2.8) 



Here is the Bose-Einstein distribution. As was discussed in [27], |28[], this problem 
can be avoided by introducing a more general and indeed more convenient definition of 
normal ordering than the one given above. It is obtained by splitting the field operators 
into two arbitrary parts \& = + (+/— do not refer to positive/negative energy 
waves) with normal ordering N ES defined to place the 'positive' parts to the right of the 
'negative' parts so that, for example 



N 



ES 



(2.9) 



where = \l/(xj), a = 1 for bosons and a = — 1 for fermions. If the contraction is 
defined in the usual way in terms of the two-point time ordered product and normal 
ordered product (|2.9|), then the operator form of Wick's theorem holds in its usual form 



provided the contraction is a c-number [Pq| . In thermal field theory, one can then choose 
a split which satisfies 

< N ES [»] > = (2.10) 

for all non-constant operators • in the interaction picture, thus mimicking the T = case 
( p. 7]) . For boson fields and all contours C used in thermal field theory |29|, this split is 
given by 



k 

53 (2^, A L)- 1 / 2 



1 



fk,n)dk,ne ik ' x + gk~A-/ kx 



fkyfldk,/!^ 



-ik-x 



+ (1 



(2.11) 
(2.12) 



where 



(1 - /fc,A)(l - 9k,ii) = 1 + iV fcj/i , fk,p9k,ii = ~ N k,ii- 



(2.13) 



Similar expressions hold for fermion field operators [p7] and also for other operators such 
as the momenta operator n. Note that N [•], the operation we call thermal normal 
ordering (TNO) when ( |2.10| ) is satisfied, reduces to : • : at T = 0. 
We will need the free thermal boson propagator 

A T {x) = < T c [}{x)4>{0)] > = e c {t) [0 + (x),0-(O)] +6 c (-t) [0 + (O),0-(x) 

where T c means contour ordering. In the imaginary time formalism (with the contour 
starting at U = and ending at —i/3), use of ( |2.11| )-( |2TT3| ) gives 



47T 



\ogjl 2 (3 2 Q 2 (x) +K+0(j2(3). 



(2.14) 



where we have expanded in powers of <C 1 and K is a constant, which will not play 
any role in our analysis (see below). The Q variable is given by Q 2 = Qq + Q\ where 



Qo{x,t) 
Qi(x,t) 



cosh(— ) sm(— , 
smh(y) cos(— , 



(2.15) 
(2.16) 



so that Q(x) is a Lorentz scalar. Observe that the boson propagator (|2.14|) is both ultra- 
violet (UV) divergent (x — > 0) as well as infra-red (IR) divergent (fl — > + ). It also 
has the same form as the zero temperature boson propagator, as for both T — > + and 
x — > + , Q a — > 7iTx ai so that the T = propagator is recovered from (|2.14j) by replacing 
Q a — > Tx a For the remainder of this section we will label the propagator also by its mass 
jl so that in (|2.14|) , A T (x) — > A T (x; ft). 

We note here other useful properties of the Q variables that will be needed later on. 



The first one concerns their short and long distance behaviour, which, from (|2.16j ) is 



Q 2 (x,i) 



(ai,t)->(0,0) 



L e 2*\*\/(> vt 



(ttT) 2 (x 2 + ?) 



and the second one is 
„ d 



with 



dx^ 
dx^ 



l —\ogZX{x) l —\ogZ Q {x) 

1 i 

logZX(x) H 



e - iei 



^ogZ Q (x) 



Zq{x) = Q {x) ± iQi(a?) 
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^log^(.) 
logQ 2 (a;) 



(2.17) 



(2.18) 



(2.19) 



2.2 The sine- Gordon and Thirring models at T > 

We now calculate the SG and MT partition functions through the following series of 
steps. 



2.2.1 Removal of UV divergences in the SG model by normal ordering 

Recall that for a scalar field theory in d dimensions with an interaction of the form r , a 
diagram with n vertices and E external lines has a UV degree of divergence D of PU| 



D = d 



E 



n 



T 

-(d 

2 y 



d 



Hence with d = 2, D = 2 — 2n so that other than the propagator (|2.14 ), the only 
divergent diagrams Vr are tadpole diagrams. Thus, as we now show, all UV divergences 
are removed by thermal normal ordering the SG Hamiltonian 



H, 



SG 



dx 



y + 2 



dx 



- — cos X(p 



7o 



H — dx 
Jo 



a 
A 
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— cos A0 + 7o 



(2.20) 



To carry out this procedure, the Hamiltonian must be divided into a free and an interact- 
ing part. Although the term cosA</> itself contains a mass term on expansion in powers 
of A, we want to keep A of arbitrary size. Consider therefore 



H 



SG 



dx [^P 2 4> 2 



d.r ( ^ cos A0 + -^P 2 4> 2 + 7o 



so that perturbations are about a scalar field of mass p. To take account of this fact, 
TNO is now denoted by iVr [•] and similarly we add an extra mass label to the propa- 
gators ( p,14| ): hence we will be dealing here with At(x; p). The link between p and the 
corresponding regularisation scale in path-integral methods is discussed in section [3[ 
TNO of ( |2.20p may be carried out by using the identity^ 



i f c d 2 xj(x)j>(x) 



e J c 



jyES e i f c d 2 xj(x)4>(x) 



e * ! c <l 2 xj(x)j)(x) 



A L d2x f c d2 yj( x ) A T{x-y;p)j(y) 



(2.21) 



Notice that had we used T = normal ordering, the zero-temperature propagator would 
have appeared in ( [2.21 ) rather than the finite temperature one: TNO means that the Q 
variables of (|2.16| ) are built in from the start. As in IJ we regulate the UV divergence of 
A T (x; p) by cutting off the theory and replace A T (x; p) by 

A T (x; p; A) = A T (x; p) - A T (x; A) 

where A is a large mass. Note that this operation cancels the additive constant K in 
the propagator ( |2.14|) . Now A^(x; p; A) is both non-singular as well as (3 independent for 
x -> 0: A T (0; p; A) = In (p 2 /A 2 ) . Also equation ( p.21| ) with j(x) = \S(x - y) gives 



A 2 



N 



ES 



so that if we define 



/ P 

a 



2 \ 87T 

A" 2 J 



then the normal ordered form of the SG potential is just as at T = [|TJ 

a 



A 



cosA</> 



— N ES 
A 2 p 



cos X(f> 



(2.22) 



(2.23) 



(2.24) 



Also TNO of Hq (which is defined in (|2.20| ) and is only equal to (|2.1| ) for zero mass 
p = 0) gives 

H = Nf[Ho\ + E T {p) (2.25) 

where, using (|2~TT| ) - (^T^) , 

MP) = U[ 



7T ' 7T 



v^l, ,2k 2 -f 

Z lil + 2N t , f) — 



2 



P 



k \ 



2k 2 



Uk„ 



4 These equalities may at first seem surprising, but they should be clarified by recalling that here 
the advanced and retarded thermal propagators are equal so that A^(x — y\ p) — At(x — y; p) — 



4> (x) , <p (y) where the positive and negative parts are given in (2.11)-(2.12). 
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and Eo(p) is infinite whereas the T-dependent part is finite. 
Thus combining (|2.24 ) with ( 2.25Q gives 



SG 



N 



ES 



Hq — — I dx cos X(j) — L"f 
X 2 Jo 



where 



7 = 70 - E T (p). 



(2.26) 



Thus ao has been multiplicatively renormalised, 70 has been renormalised by an infinite 
temperature independent part appearing in Et(p) as well as having been shifted by a 
finite temperature dependent amount, and A is unchanged. 

Thermal normal ordering the SG Hamiltonian has therefore absorbed all UV infinities 
just as zero temperature normal ordering does but it has also introduced some extra 
T-dependent finite terms. Its power will become apparent in the next section when 
thermal expectation values are calculated. 



2.2.2 Perturbation theory in the SG model 

We now remove the IR divergences of the thermal boson propagator by introducing a 
mass ft, as in Q2.14| ), into the SG Hamiltonian (|2.20 ). At the end of the calculation we 
take ji — > + and hence are free to add the extra mass term within the normal ordering 
giving the Hamiltonian 



H 



N ES 
P 



N ES 



H + / dx 



1 ~2 2.2 a 



- — cos X4> 



(A§ - 7Z,) - J dx (cos A0 



A 7 



ES 



a 



dx (cos 



= A + Aj 



Here A and Aj denote the free (a = 0) and interacting Hamiltonians respectively. 

As usual in perturbation theory, the thermal expectation value of Heisenberg operators 
•h are written in terms of interaction picture operators • as 



< *H > 



Tr{< 



-PH. 



Tr {eH 3 ^} 



Tr [e~P Ao T c 


u{u- 


-iP,U>]} 


Tr {e-^»T c 


U(ti 


-0,U)]} 



with times lying on the contour C. For = 

U(-i/3,0)=T c 



i f ^dt"Aj{if') 



Thus on expansion in Aj and working in Euclidean space, the SG partition function is 
given by 

Z SG (T) = lim V 1 f dh . . . f dt n Tr {e^ A °T c lA^h) . . . Aj(t n )} } , (2.27) 
nl Jo Jo L L J J 



10 



and one is left to calculate free expectation values; 



Tr | e -Mo # | 



(2.28) 



Here the prime label occurs because the operator appearing in the thermal weight is 



N 



ES 



n 



and not the simple diagonal operator : Hq Y^k u k,^\ fflKfi as m (EH)- 
Since the results of p7|| can only be used with this latter thermal weight, we need to 
manipulate <C^o- For that purpose, note that 



ES 



so that from (12. 31) one has 



A 7 



ES 



n 



2M 



fe, P j 



' 2F + p_ 2 + A 



Pi2 > 



Hence substituting into ( |2.28| ) gives 



< • > 



and the results of |27| may then be used to say that <C N? s [•] ^>o= whenever • is a 
product of operators (not including a constant). Finally, before proceeding, we note from 
equation (|2.22| ) that 



N 



ES 



(2.29) 



Following the steps in [ffl] for T = 0, let us calculate at T > the free correlator 



« r c n at 



ES 

p 



>r 



Tr < 


' e -f3N? S [H§] Tc 


11'/ : \ 




]} 


Tr< 







Tr {< • "• : -V r [n™ =1 N? s 



where we have used ( |2.29| ). Note that the above correlator is the one we would need to 
prove the equivalence ( |1.3|) by expanding in a 0. Also from Q2.21|) with j(x) = Xj5(x—Xj) 
and using the fact that time ordering always puts the operators in a given order which 
depends only on their time arguments, one has 



ES 



e §A T (0;A)£,.A 



2 

3 ' 'j 
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Hence 



3=1 



ES 



N ES 



N ES 



iV, 



£5 



g^EjOife)) 



g-l E, E fc ^Ar(*j-**;A)A* e | E, ^ 2 At(0;A) 



"Ej>fc ^]&T{Xj-X k \fl)\ k 



nn[^i-**) 

j=lj>k 



A j A fc 



giving, finally, 



3=1 



>r 



V A 2 



ss 1 



* n n [^wo^ - **) 



>o x 



' rj2 \ 87r ^ )? A j 71 



|2l 4ir/ 



nn[^ 2 iQ(^--x fc )i 2 j -(2.30) 

\V J j=lj>k 

where the last line follows from expanding the exponential inside the normal ordered term 
and noting that all terms which contain cf) n vanish (by fl2.10|) ) apart from the one for which 
n = 0, giving a contribution of 1. It is at this stage that thermal normal ordering has 
been crucial for the calculation — had we used zero temperature normal ordering, this 
correlator would have been very much more difficult to calculate. Note that the terms 
proportional to ft in ( |2.30|) have a contribution 



V a 2 n \ \ 

3=1 i>fc 



(En* 
\1 4lr 



(2.31) 



so unless En \ = 0, the result vanishes as p, — > 0, the limit in which we are interested. 
Hence we only consider configurations with J2 n <\? = 0, a condition which will become 
analogous to the fermion chiral selection rule of section |3|. In ( |2.30|) n must therefore be 
even. If we let m = n/2, yj = Xj for j = n/2 + 1, . . . , n and define A± = N ES 
then 



< T c 11 A+ixjjA-iyj) > = 11 - 2 

nr=i[/?viQ(^--^)i 2 ]- 

To calculate the partition function we need 



2n 



3=1 



1 (2n)! " nj >fc [/?V |Q(*, - x k )\ 2 \Q(y 3 - y k )\ 2 ]^ 

2 11 



2 2n 



A 2 

2 S 



(2.32) 



which follows from (|2.32|) and ( p.31|) . Observe that the total number of f3p factors in this 
expression is (f3p)~^( 2n \ Hence the partition function (|2.27|) is 

Z SG (T) =Tr{e-^ G } = 

12 



2/i 



77! 




d Xi . . . d X2n <C 2c 

2T! 



2h 



IJ 5 [cos A0(tj 
7=1 



>0 



n 



U] >k [\Q(x,-x k )\ 2 \Q(y 3 -y k ] 



A 2 

I215F 



I 21 4vr 



(2.33) 



Tr{< 



-Mo | 



e 13101 Zq B (T) and it is finite for 70 satisfying (|2.26|) with 7 



where Z S (T) 
finite. 

One might wonder whether there are any extra divergences in the Xj, yj integrals 
in ( |2.33|) , which, for a given n correspond to closed diagrams with In vertices and an 
arbitrary number of loops. Recalling the asymptotic behaviour of the Q variables in 
one realises that the integrals in (|2.33|) are always finite for large spatial Xj or yj, 
because the number of Q 2 (xj) powers in the denominator is always bigger than that in the 
numerator. As for the behaviour of the integrand when Xj — > we see that the integral 
converges provided A 2 < Att. Hence in this case, all the UV divergences of the theory 
are removed by renormalising a. If A 2 > 4tt, new divergences could appear in principle. 
However, it can be shown that as long as A 2 < 8tc, those extra divergences can also be 
removed, by subtractions to the free energy |5I |. The same argument holds in our case, 
since the UV divergences are T independent. In fact, for A 2 < 8tt the theory is super- 
renormalisable [|] . For A 2 > 8n though, the theory is no longer super- renormalisable and 
indeed the ground-state disappears from the spectrum [|TJ. These bounds should borne in 
mind and, although formally our expressions and results hold for any A, they will remain 
finite only for A 2 < 8tt where we will assume that the T = subtractions have been 
already made. 



2.2.3 Perturbation theory in the massive Thirring Model 

In the Thirring model, the situation is somewhat unusual as perturbation theory not 
about the zero coupling g = theory, but about the massless itlq = theory. The 
correlator corresponding to (|2.32| ) is therefore 



«T c n^(^)^(%)» 
3=1 



where R stands for renormalisation and 



m = 
9 # 



(2.34) 



a±{x) = tp(x)P±ip(x) with P± 



^±75) 



At T = when the thermal expectation values are replaced by vacuum expectation values, 
such correlators have been calculated in a canonical operator approach by Klaiber |3"2"f . 
His techniques have been used extensively by Coleman [p]] as well as by others for example 
in the study of the Schwinger model [0. However, the generalisation of Klaiber's work 
to T > does not seem to have been carried out. As in the case of the SG model, 
this is obtained by firstly replacing zero temperature normal ordering by thermal normal 
ordering and then using thermal expectation values rather than the vacuum expectation 
This changes all zero temperature propagators to their finite temperature 
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values of 32 



analogues just as in section |2.1.2| , and is equivalent to replacing x a —> (3Q a . Furthermore, 
the regularisation remains unchanged from the T = case — that is, given the parameter 
5 defined in [|32]|; 

5 _g 2 2* + g 2 

477 7T + g 2 ' 



it follows from equations (VII. 5), (VII. 4) and (VII. 1) of [32] that in order for the correlator 
of <j±(x)'s to be finite, these operators must be regularised as 



(2.35) 



Here the ^j's label the components of the spinor, e' is a 2- vector, and M in is an arbitrary 
renormalisation scale which, as opposed to has been included so as not to alter the 
dimension of the operators. Finally, <r^ is given by the hermitian conjugate of ( p.35|) . 



As we will verify also in section |3.3| , following these steps and using the equations of |32 
mentioned above, one obtains 



< T c J] a^x^iyj) > m = 
i=i 



9 ^ 



M 2n U n [^M 4 |Q(x, - x k )\ 2 \Q( yj - y k )\ 2 ]^ 

(?) 5 = — 



UU [/3 2 M 2 |Q(x,-rf 



(2.36) 



Observe that up to the factor of (M/2) 2n , equations ( |2.36|) and (|2.32|) are identical at 
T > as long as the two theories are renormalised at the same scale, M = p, and provided 
that 

A 2 



7T 



An 



which is just ( |1.5| ). The partition function then follows directly, and is given by 

OO 1 2 T ^ / "I TL 



Uk<j [Q 2 (xj ~ x k )Q 2 (y j - y k )} 



l-K 2 /n 



UU [Q 2 (^-y k )] 



1-k 2 /tt 



Here Zq is given in 



l + g 2 /7r 

and m is the renormalised mass (see also section 0) 



m = m (s'p) 



-26 



(2.37) 



(2.38) 



(2.39) 



In perturbation theory the SG and MT model partition functions are therefore identical 
term by term if we identify 

a 

pm = — 



as in (|1.6|). 

In order to prove the relationship between arbitrary products of zero charge operators 
at different space-time points in the two theories and treat the /i / problem, we will 
turn to more powerful path integral methods. 
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3 Bosonisation with Path-Integral methods at T > 

and ji = 

Again we work in path-integral imaginary time formalism of thermal field theory and 
hence in Euclidean space-time with metric (+,+). Now we change notation slightly so 
that t = i or x° denote the (real) Euclidean time, but, unless stated otherwise, the 
remaining notation follows that of pl| . As usual the exponential factor in the path 
integral is exp(— S) = exp(— J T cPx£(x, t)) with J T d 2 x = Jjf dt f^™ dx, and the Euclidean 
Lagrangian densities for the SG and massive Thirring models are respectively 



-SG 



\ 9 * 



- — cos X(p, 
1 



-V>(0 + m )ip + ^g 2 j^(x)f(x) 



where we set 70 = in the SG Lagrangian. The Euclidean gamma matrices are 



7 




7 







75 



so that 



V7 1 




with eoi = 1. 

Let us briefly sketch the outline of this section. Firstly, in section 3.1 , some useful 
free boson and fermion correlators are calculated in the path integral formalism for T > 
and fi = 0. Here we make use of the free generating functionals and thereby reproduce 
the results for those correlators obtained with canonical methods in the previous section. 



The aim of section [3.2| is to recall how renormalisation takes place in the path integral 
and, in turn, to compare with our renormalisation conventions in the canonical formalism. 



Finally, in section [O] we re-obtain the equivalence between the partition functions of the 
two models, now from the path integral viewpoint. This then serves to introduce the 
SG and MT generating functionals which will be needed later to re-derive the correlator 
( p.36| ) in this language, as well as to introduce the auxiliary field technique we will be 
using extensively. Besides, most of the path integral techniques and results needed in the 
rest of the paper will also be introduced in section |3.3| . As in the previous section, we 
always proceed by expanding around the ckq = m o = theories. 



3.1 Free generating functionals and useful correlators 

We start from the free boson Euclidean generating functional 



Z^[J;T] = 







/ d(j) exp < — 


[ d 2 x- 


'periodic I 


Jt 2 L 



(<9 Q 0) 2 + /lV +J(x)(f>(x) 



Z B [0;T]expjijT (fx jT d 2 y J(x)A T (x - y) J(y)} (3.1) 



where the propagator is given in (|2.14 ) for small /it/3 and the free boson partition function 
is Zq(T) = Zq[0;T] as in ( |2.4j ). Note that from this section onwards the p, labels in 
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propagators and partition functions will be removed since there is no possible confusion 
with any other mass scale. Similarly the free massless fermion generating functional is 



Z F [fj, V ;T] 



Nq I dip dip exp 

J antiper 



d 2 x 



+ ijip + ipf) 



Z F [0,0] exp 



d x id yr](x)S(x — y)r](y) 

T JT 



(3.2) 



where Z F [0,0] = Z F in ( |2.5|) and N F is the fermionic counterpart of Np in (|1.7| ). The 
free massless propagator S(x — y) can also be expressed in terms of the Q variables as 



24 



1 Qopiz) 
2/3 Q 2 {x) 



(3.3) 



S a(3 (x) -- 

where the indices a, f3 are Dirac indices. 

The path-integral counterpart of the free correlator ([2.32 ) is easily calculated within 
the path-integral formalism. We outline the main steps as they are typical of those 
required below. We define now A± = exp [±iA</>] (remember that there is no operator 
normal ordering in the path integral), so that, 

7B 



11 7 B \ J- Tl 

« T c J] A+MMVj) »o= -fw?^ 
j=l z o l J ) 



with 



J{z) 



n 



2), 



6®i 



z - Vi> 



3=1 



Unlike section ^, here the UV divergence of the propagator ( [2.14 ) will be regulated by 
replacing 

Q 2 (0, 0) -> Q 2 (e , £l ) = T 2 e 2 + 0(e 3 ) (3.4) 



where e a — > + and 



7T 



(el + ef). Then from Q, (|1|) and O it follows that 



nLi [Q 2 (xj-yk)} 



Observe that this correlator is ill-defined as e — ► + . Since we are still dealing with 
the free theory, the parameter A entering the definition of the A± operators above is, so 
far, arbitrary and, hence, the above correlator could, for instance, diverge or vanish by 
taking A purely imaginary or real respectively. This behaviour is due to the short-distance 
(UV) divergent behaviour of the composite operator exp[i\(j)(x)} which needs then to be 
renormalised (see section |3.2.1| below). Also note that the ft, dependence in ( ^.5[ ) has 
vanished for the same reasons as those of section [2.2.2| , i.e, because we have chosen equal 
numbers of A + and A-. It is not difficult to see that, had we chosen different numbers, 
the above correlator would have been proportional to ft a with a > and hence it would 
vanish in the ft — ► + limit. This is the bosonic counterpart of the so called fermion chiral 
selection rule which we will discuss below. 

Similarly we need to calculate the fermion free correlator of an arbitrary number of 
insertions of the operators a±(x) = ip(x)P±ip(x). To do that, recall that the massless 
fermion theory is invariant under chiral transformations ip — > explia^ip with a real 
arbitrary. Under such transformation cr±(x) — > exp(±2ia)a±(x) and therefore the ther- 
mal average of a product of <J±(x) operators will vanish in the massless case unless the 
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number of <7+ and <r_ is the same. This is the chiral selection rule, which only holds for 
mo = 0. Following [21], the required correlator is obtained by shifting if) — > t/>7° so that, 
naming i\) a with a = 1, 2 the two components of the bispinor i/j, the free massless theory 
now decouples into two free theories for the spinors ip a , with propagators 



S a (x - y) = /r 1 J2 [ 

n&TL 



dpi e _ ip . (x _ y) ipo±p 



1 



1 



2tt 



2(3Z±(x-yY 



(3.6) 



The +, — signs correspond to a = 1, 2 respectively and Zq are given in ( |2.19 ). Now the 
transformation ip — > -^7° takes a + — > ^2^1 and er_ — > ip\ip2- Recall also the thermal 



averaged version of Wick's theorem for fermions |21[ which is 



< T c Y[ i) a {xj)i)a{yj) >o= det S a 

3=1 



(3.7) 



where the subscript for fermion correlators will denote the free massless case m = g = 
and S a denotes the n x n matrix with elements (S a )jk = S a (xj — i/k)- Hence, using the 



result 21, 24 



\n+l 



det 



Z Q {xj-yk) 



n 

3=1 



Uk<j Zq(xj - x k )Z±( yj - y k ) 



I 0C 4 



Vh) 



(3.8) 



where Zq is given in ( |2.19| ) we have 



« T c J] a+ixja-fa) » = (2/3)' 2n n 



IL<j [Q 2 (^j - x k )Q 2 {y j - 



(3.9) 



nLi[Q 2 (^-^)l 

Notice that the above correlator has exactly the same structure as the boson correlator 



(p.5|). However, (|3.9|) is finite and well-defined, unlike (|3.5|) , since it contains no products 
of fields at the same point and also because there are no mixing terms between ip\ and 



ip2 in the Lagrangian. We will come back to this point in section 3.2.2 



3.2 Renormalisation 

Here we recall the path integral renormalisation procedure for our present case. We want 
to remark that our regularisation and renormalisation conventions are going to be the 
same as those discussed, for instance, in ETJ for T — \i — 0. That renormalisation is 



independent of T and fi is an expected result which will be explicitly checked throughout 
the rest of the paper. In addition, note that we need only to analyse renormalisation in 
the free («o = 0) SG model and in the massless (m = 0) Thirring model. The reason 
is that, just as with the canonical methods, we treat the SG model and the massive 
Thirring model always in terms of free and mo = correlators, with an arbitrary number 
of insertions of (ao/A 2 ) cos cf>(x) and moTpip^x) respectively. 



3.2.1 sine-Gordon model 

As commented above, the ill-defined behaviour of the free correlator (|3.5| ) is caused by 
the short-distance behaviour of the free propagator in (|3.4j) . This problem can be solved 
by renormalising the composite operator exp[ia<j)(x)] with a arbitrary as 



[exp [za0(x)]] 



bare = [«P M(x)]] 
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R 



(3.10) 



where p is an arbitrary renormalisation scale. Note that this expression is analogous to 
( p. 22 ) with the identification 

1 

£ = v 

although in the operator formalism the renormalisation was carried out through TNO. 
Also observe that the dependence with the scale p is the same, though it has slightly 
different origins — in the canonical operator approach p corresponded to an arbitrary 
mass at which normal ordering was performed whereas here there is no normal ordering 
and p is an arbitrary renormalisation scale. 

As observed in section [2.2.2| , all the infinities that appear in the theory are just those 



coming from ( |3.10| ), after free-energy subtractions, provided A 2 < 87r. Those divergences 
can be removed by defining the renormalised coupling constant a as 

a = a (ep) x2/4n (3.11) 

(c.f. equation ( |2.23| )), and the coupling constant A does not need renormalisation. We 
will check explicitly below that with this coupling constant renormalisation all our results 
remain finite (for A 2 < 8n, as explained above). 

3.2.2 Thirring model 

We have seen that the u±{x) correlators are finite in the m = g — case because the two 
components of the bispinor do not mix in the Lagrangian. If g ^ 0, m = 0, the chiral 
symmetry is still unbroken so that u±{x) correlators again appear in the same number. 
However, when ip — > ^7°, (^jf^if)) 2 — > / ^4>\'^>\^)2'4 } 2 and therefore there is mixing between 
■01 and ip2- Hence 

< M x )Mx)Mv)Mv) >^ - < M x )Mv) >< Mv)M x ) > 

when g ^ 0, so that products of fields at the same point appear and then the u±{x) 
correlator becomes divergent. As for the boson operator exp[ia<f)(x)], the <J±(x) composite 
operators need renormalisation. Also as in SG, those are the only infinities to consider 
and now they can be absorbed in the renormalised mass m. We anticipate here the results 
for such renormalisations, which we will explicitly check below, 

[a ± (x)} R = M" 2/ >±(*)] 6are (3.12) 
m = m (ep)- K2/n (3.13) 

where both models are renormalised at the same scale p and k is given in ( j2.38j ). Notice 
that our renormalisation conventions here for the fermion case differ from those in section 
2.2.3| , which are those used in (compare for instance (|3.13|) with (|2.39|) ). For that 



reason, we have used different notations for our short-distance cut-offs, namely, e and e'. 
However, observe that from ( |3.13|) and (|2.39|) and the definitions of K and 5, they are 
related through 

e = e (ep) 27r +<>' 2 

and since the exponent of e in the r.h.s of the above equation is positive for A 2 < 87r 
according to (1.5|), we can take both e,e' — > + consistently. In addition, we will see 



below that the answer for the (finite) renormalised correlators and partition functions in 
both canonical and path integral approaches is the same, even though our renormalisation 
conventions are different, which is a necessary consistency check. 
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3.3 Generating functionals and partition functions 
3.3.1 The SG generating functional 

As usual, all correlators can be obtained from the generating functional which, in the SG 
model, is 

Z SG [J; T) = Np ! d<P exp {- [ / d 2 xC SG [<P} + J{x)<f>{x)\ } (3.14) 

J periodic I [jT J J 

where Zsg{T) = -Zsg[0;T] is the SG partition function. Expanding formally in powers 
of cuo / A 2 gives 

ZscU;T] = E ^ E ( I ) U.J T d^x 3 z»[J + J;T] 



with 

n 

J( z ) = -i\^ Sj 5W( z - Xj ) 

3=1 

Hence, from (13.11) we have 



a 



+ j = l,...,k 

j = k + 1, ...,n 



Z SG [J;T] = Z^[J;T]Y: 



n=0 



n\ V2A 2 



E ( J n/ T ^Sexp|-^nA T (0) 

n „ n I 

— iA E s j / d 2 xAT(xj — x) J(x) — A 2 E E s j s kA>T(xj — Xfc) > . (3.15) 
i=i 7t j=ifc<j " J 

The above expression still contains both the UV (e — > + ) and IR (Jl — > + ) divergences 
of the propagator. However, we will see that the observables and correlators of interest 
are finite, after renormalisation, when those limits are taken and furthermore that the 
p, —>■ + limit selects a particular value of k in ( [3. 15 ), which (see below) is again equivalent 
to the fermion chiral selection rule. 



3.3.2 The MT generating functional 

The MT model generating functional is defined, as customarily, as 



Z M1 [fj, 77; T, (i]=Ng / exp 

J antiper 



(3.16) 



where, for completeness, we have also included a chemical potential /i term that we will 
analyse in section |5|. Following || |3], [21|, |24|, let us rewrite the quartic Thirring interaction 
as a "gauge-like" interaction with an auxiliary vector field 



z MT [fj,v;T,n] 

with 

Zf[A,fj,ri',T,n] 

and 



Z /\ ( I ) J periodic 



dipdip exp 

antiper 



[ dA^e W \~ [ d 2 xA^x)A"(x)}z f [A,fj,r];T,fx} (3.17) 

J periodic V Z JT 



d 2 x Up ^ + i<7 4( x ) + m o — ^7°) if) + ryij) + tpry 



(3.18) 



Za{T) 



dA^ exp 

periodic 



2 7t 



d 2 xA ll (x)A^(x) 



In the remaining of this section and in section [| we will concentrate only in the case 
/i = 0. The /i 7^ case will be studied in section || 
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3.3.3 The partition functions 



We now turn to reproduce the result in |24|] for the partition functions of the SG and MT 
models, paying particular attention to UV renormalisation and to the \x — + + limit. 

Let us start with the SG partition function, which is obtained by setting J = 
in Q3.15I) . Using the explicit form of the propagator in (|2.14 ) one finds that Z$g is 
proportional to jl a with 

\ 2 n \ 2 

« = -(E^) 2 = 7r(2fc-^) 2 



An 



An 



Hence, in the jx — > + limit the only term that survives in the fc-sum is k = n/2, which in 
turn selects only even values for n in the n-sum. Here we also see why, in the path-integral 
language, the value of K in ( |2.14f ) is irrelevant: the reason is that any constant term in 
A T cancels when J = in (|3.15|) since it is multiplied by (2k — n) in the only surviving 
term k = n/2. Then, redefining variables for j = n/2 + 1, n and n — > 2n 

and taking into account the definition of the renormalised coupling constant a in ( p. 11 
we get 



Zsg(T) 



n=0 



1 ' 

n\, 



i 


a 






2>? 





2/i 



X 



n 



d 2 Xjd 2 ijj : 



U k <j [Q 2 (xj - Xk)Q 2 (yj - yk)] 



A 2 /4tt 



YIU [Q 2 (xj-yk)) 



(3.19) 



(c.f. equation ( |2.33| )). As in section p.2.2| , we again see that, at least for A < An, all 
the UV divergences of the theory are removed by renormalising a (see our comments in 
section |2.2.2] about the regime An < A 2 < 8n) . 

Our next step will be to calculate the MT model partition function and compare it 
with ( |3.19| ). First, we write the auxiliary vector field A^x) as 



21 



A^(x) = — [<V?(x) + ie^ v d v (f)(x)] 
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(3.20) 



with r)(x) and 4>(x) two boson fields, chosen to be periodic in Euclidean time. Then, by 
performing the following chiral transformation in the path integral fermion variables 

ip(x) = exp {i [n(x) + 750(x)]} ip'(x) 
4>(x) = ip'(x) exp {i [—n(x] + 7 5 0(a;)]} 



the /i = MT partition function, obtained by setting n 
may be rewritten as 



v = n = o in (gl7D - (PH) , 



J MT 



(T) 



Za(T^ J antiper J periodic 



+ 



2^ [ ^ 



dipdi/j / dcf) dn exp < / d 2 x 



2iy 5 if>(x) 



1 



(x)] 2 + ip(x) $ + moe 



2g 2 

ip(x 



(3.21) 

Here J a is the Jacobian of the transformation dA^ = J^d^dr], k is given by ( |2.38| ) and we 
have used that the Jacobian of the change of variable dipdip = Jpdip'dip' , which is given 
by the axial anomaly 



33, 21 



J* 



exp 



7^ / d 2 xd^(x)d^(x) 
In Jt 
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Finally we have relabelled the fermion variables ip' — > ip. Note that in ( |3.21 ) the 77 field is 
free and the only interaction term is the one proportional to mo which mixes the fermion 
fields and <f>(x). Now expand in powers of mo, note that ipe 2l ' y5 ^ip = a+e 21 ^ + er_e~ 2 ^ and 
recall that, due to the fermion chiral selection rule, only terms with the same number of 
a + and cr_ contribute. The integral over rj is readily done and the boson and fermion 
correlators appearing in an arbitrary term in the expansion are of the form ( |3.5| ) and 
(p.9|) , though when using ( |3.5|) in (|3.21|) the fields must be rescaled — > intfi because of 
the (1/k 2 ) factor in the kinetic term. Finally, recalling that Za(T) = Ja(Zq(T)/Np) 2 up 
to T-independent terms and then relabelling the sum and integral variables gives 

2n 



J MT 



(T) 



n=0 



m 
2P 



K / TV 



X 



n 



d 2 Xjd 2 yj ^ k< ^ 



[Q 2 {x j -x k )Q 2 {y j -y k )\ 

1 1 — K 2 /7T 



1-k 2 /tt 



(3.22) 



nju [Q 2 {^-vk)Y 

where we have renormalised m as ( |3.13|) . This is just equation ([2.371) an d shows how the 



work in simplifies in the path integral framework. Again by comparing fl3.22| ) with 



( p,19| ) and using ( |2.6|) we reach the same conclusion as that of section [2.2.3| and 

Z MT {T) = Z SG (T) (3.23) 

up to vacuum energy and T-independent terms, provided the renormalised coupling con- 
stants of both models are identified as in ( |1.5|) -( |L6[) . 

At this point, notice that we could have expanded directly in mo in the original MT 
Lagrangian, and then every term in the expansion would have been the g 7^ 0, m = 
correlator of a product of <j±(x). Thus from (|3.22j ) one obtains 

"(re) -' 



«r c n<7+(i>-y > 



m = 
9 ^ 



- K 2 /tt 



2(3 

A Uk<j [Q 2 (xj - x k )Q 2 ( yj - y k )} 1 ~ K2/7 ' 



which, as we said in section [3.2.2| , is infinite. However, with the renormalisation Q3.12 ) 
for the composite operators <J±(x) one obtains ( p.36| ) with M = p and so this is just the 
T > extension of the result in 0] , as required (see our comments in section |3.2| about 
our renormalisation conventions). 

Also once again comparison of (|3.24 ) with ( |3.5|) gives for the renormalised correlators 



3> m = 
9^0 



2ll 



«r c n^(^(%)» 
3=1 



(3.25) 

(p/2) exp(±iA0). This is the 



as long as (|L5| ) holds, which suggests the equivalence a± 

m = a = version of (|1.3|) , just by combining it with ( |l.6|) . In the next section we will 
obtain the m^O version of this formula. 



4 The mass and current equivalences at T > and 

li = 



In this section we show that the mass and current operator equivalences (1.3)- (|1.4j) also 
hold at finite temperature, i.e, for thermal averages of correlators of those operators. The 
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condensates of a± and A± are analysed in section |4.1| whereas we consider an arbitrary 
chiral invariant chain of those operators in section i2. The thermal current equivalence 
( |1.4| ) is discussed in section [D| where we make use of point-splitting regularisation. 



4.1 The mass equivalence for the condensates. 



The simplest correlators to look at to prove would be <C A±(x) ^>sg and 

<C cf±{x) ^>mt which do not vanish for m ^ because the chiral symmetry is bro- 
ken. However, their relationship is actually a consequence of the partition functions 
equivalence ( |3.23| ). The reason is that they are independent of x simply because of the 
translation invariance of the theory so that 

< [cosA0(a;)] 6are > SG = < [cos \<j>{x)] R > 5G = Jl a Q^ Z SG{ T ) 



m -C 



bare 



> MT = m -C 



tfjlfj 



1 d 

_ >a/t = —m—Z MT (T). 
R pL dm 



In fact, the above fermion correlator is nothing but the fermion condensate which plays 
the role of the order parameter for the chiral phase transition. However, we do not expect 
chiral symmetry restoration at finite temperature here since Coleman's theorem prevents 
any continuous symmetry to be spontaneously broken in 2D |20|| . 

Therefore to prove the desired equivalence for the above correlators, all we have to do 
is to differentiate in ( |3.19| ) and (|3.22| ) with respect to a and m respectively. However, it 
is instructive for the next sections to carry out an explicit calculation. From (|1.7| ) and 
( P-14| ), we have 



< A± > 



SG 



z SG [J±;T} 

Z S g(T) 



with J±(x) = =rA5 (2) (x - z). 



Next use ( |3.15| ) to obtain the selection rule for this correlator: the total power of ft is 
now given by (A 2 /47r)(l ± 2k — n) 2 , so that only k — (n =p l)/2 and hence odd values of 
n in the sum, contribute now to (|3.15| ). Then, by redefining n — ► 2n + 1, we obtain 



< A±(x) > SG 



x 



Zj(T) 
Zsg(T) 



E 



n\(n + l)\ 



n=0 

n n+1 



a_ 




2A 2 


\P / 



2n+l 



x (Te) A2 ^nn lj\d 2 Vl 

3=1 1=1 



T 



x) 



Q 2 (yi - x) 



A 2 /4tt 



n fc <j U k '<i [Q 2 (xj - x k )Q 2 (yi - yk')] 



[Q 2 (^-yi)] 



A 2 /4tt 



(4.1) 



Here we have relabelled Xj *- 



yi and used Q 2 



[x 



Q 2 



-xj 



Observe that A±(x) ^> is the same for both ±, which we could have expected by 
realising that both the SG Lagrangian and the path integral measure are invariant under 
4> — > —cf). Also the UV e dependence in (|4.1|) can be absorbed in the renormalisation of 
exp(iA0) according to ( |3.10|) so that once again there are no more infinities for A 2 < 4ir 
(see comments above). Thus we are allowed to shift the integration variables as Xj — ► 
Xj + x and jjj — > y d - + x Vj, which shows explicitly that the above correlator is actually x 
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independent. The final answer is 

Z*(T) 



E 



a 



2n+l 



Z SG (T)^ n!(n + l)! V2A 2 , 

(n+l)A 2 /27r n n+1 

n n / d2x ^ 2 yi 



3=1 1=1 



Q\x 3 



-i A 2 /4tt 



Q 2 (yj) 



X 



IIfc<j n^<j [Q 2 (xj - x k )Q 2 (yi - yk')} 



A 2 /4tt 



(4.2) 



which vanishes in the free case a = 0, consistently with the discussion in sections |2.2.2 
and CO. It can also be checked, order by order in the a expansion that the above formula 
indeed gives the same answer as (\ 2 //3L)(dZsG/da), as explained above, which is a 
consistency check. 

Now calculate the a± (x) correlator in the massive Thirring model. Following the same 
steps as those leading from ([3.171) to (|3.21|) and taking into account that 



ip exp[2i7 5 0]P ± ^ = a±(x) exp[±2i(j 



gives 



< G ± (x) >MT 



Z F (T) 



Zmt(T) ^ n\ 



oo n n 

E^fE 



n 
k 



k n—k 



k=0 
k n—k 

x < T c a±(x) J] J] a + ( Xj )a4 yi ) > 



t\i it n, n 

WW Ux jd \ 

3=1 1=1 J1 



j=l 1=1 

k n—k 
j=l 1=1 



>0 • 



The chiral selection rule then forces the choice k = (n=p l)/2 and n > 1 odd so as to have 
the same number of cr+ and cr_. Now let n — > 2n + 1, use (3.5) and (3.S) and renormalise 
the mass and the composite operator <J±(x) as in (|3.12| )- (|3.13|) . One obtains 



Zjf(T) 



1 p / mp\ 
Z MT (T)^ n\(n + l)\2 \Y ) 



oo 

E 



2ra+l 



-2(l-K 2 /7r)(n+l) 



n n+1 „ 
7=1 2=1 J7 L 



X 



X 



1-k 2 /tt 



.Q 2 (^)_ 

n fc <j E[fc </ [Q 2 {xj - x k )Q 2 {yi - yk')} 1 '^ _ 



1-k 2 /tt 



[Q'^-yi)} 

Comparison of Q4.3|) with (|4.2|) and use of the equivalences (|l~5| ), (|L^ 
finally gives 

< af (x) > A fT = f < ^±(&) >5G • 



(4.3) 
D and (p|) 



(4.4) 



As this identity depends explicitly on the scale p, it can be combined with (Of) to get 



- - -I a 

i)(x)i)(x)\ r > MT = — < [cOsA0(x)] fl >5G 



(4.5) 



which is scale independent. 



23 



4.2 The mass equivalence for chiral invariant operators 



We now prove ( 1-3 ) for thermal correlators of an arbitrary number of a± operators eval- 
uated at different space-time points. As was commented previously, there is no need 
to consider the same number of a + and cr_ in the m 7^ case. However, such a chiral 
invariant combination is chosen here because then the selection rule will be just the same 
as for mo = which clearly simplifies matters. We do not present the calculation for 
more complicated cases here since it is more involved but it goes through in the same 
way and leads to the expected equivalence, as we will comment later. 
Once again we start from the SG model and calculate the correlator 



N 



A N (X; Y) =< T c J] A+iXAA-M) » 5G = ^ff^} 
3=1 AsgU J 



with 



N 

j N (z) = -i\j2[s {2) ( 



z - X*\ 



(4.6) 



Notice that differentation of the partition function would only give global information 
about the above correlator, i.e, its total integral over Xj and Yj, which is not enough 
to prove the required equivalence. As before, on substitution of ( |4.6| ) into ( |3.15| ) only 
k = n/2 survives. On the other hand, Zq[J n ,T] gives precisely the free contribution in 
( ft^D so that 



A R N {X;Y) 



§^A» N (x-,Y )o ±(r 2 



n=0 



nl 



_a_ r 

2A 2 L 



A 2 /4tt' 
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A 2 /4tt 
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X 



v r Q 2 (x j -X k )Q 2 (y j -Y k ) 
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k=l 



Q^ Xj -Y k )Q^ yj -X k 



A 2 /4tt 



where, according to ( |3.5|) and (|3.10| ). 



^{X-Y\ 
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\ N\ 2 /2tt N 



n 

3=1 



n fc<J [Q 2 (A,-x fc )Q 2 (y,-n)] 



A 2 /4tt 



nti [Q 2 {x 3 -Y k )} 



A 2 /4vr 



(4.7) 



Now, in the Thirring model we have 
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« T c J] a^X^iY,) »mt 
i=i 

r 7F(rp\ OO 

M 2 ^#^gr E 



00 2n 
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(T) 



ni 



n+N 

n 

i=i 



d 2 Xjd 2 yj 
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Y[5^(x k -X k )6^(y k -Y k 



Lk=l 



n=0 

n+N n+N 

x « T c 1] a + (x 3 )a_(y 3 ) » « T c JJ e' 2 ^^^ 

3=1 3=1 
yF(rp\ N 00 , I s 2 

= zm <<c Tc n ^>-(^) » »=; E Q) 
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k 2 /tv' 



2n 
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X 



II d 2 x ] d 2 y j 

3=1 T 



Uk<j [Q 2 (xj - Xk)Q 2 (yj - yk)] 1 k2/7T 



UU [Q 2 (xj-y k )] 



1-k 2 /tt 
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n 



1-k 2 /tt 



Q 2 {x 3 -X k )Q 2 (y J -Y k ) 

k AlQ 2 ( x J- Y k)Q 2 (yj-Xk) 

where, from (gj^) and ( gig) 

iV 

«T c n^>^)» 



(4.? 



m = 
9 ?S 





( py 2 /*' 


2/3 





2N 



x 11 — 



— • (4.9) 

Therefore, using ( |3.25|) and the equivalences between the coupling constants and partition 
functions, we find the desired result 



« T c II a+(Xj)(T R (Yj) » MT = (g « T c J] ^(X,)^(r,) » 5G . 
j=i w i=i 



For more complicated strings of cr± the calculation goes through in a similar way and 
leads to the expected result in which o±(xj) — > (p/2)A±(xj). 

4.3 The current equivalence 

We now analyse the current equivalence (|1.4j) . When dealing with correlators including 
insertions of the current operator j^(x) in the Thirring model, there is an added compli- 
cation since a product of fields at the same point gives a divergence. We thus have to 
specify how to regularise and, eventually, renormalise. To do that we follow the same ap- 
proach as in |19|, |25[ for the massless case and use point-splitting regular isat ion. Indeed, 



as we will expand in the fermion mass, it will be enough to regularise the current for 
m — 0. Other than that, the previous discussion about the renormalisation of the theory 
will remain unaltered, i.e, no further infinities will appear once we choose the proper 
regulator (see below). We begin by quickly sketching how the point-splitting method is 
applied to this case. 



4.3.1 Point-Splitting regularisation for the fermion current 



The point-splitting regularisation method 34 consists in replacing products of fermion 



fields ip a (x)ip/3(x) — ► \im s y ^ x ip a (y)ipp(x) where lim*^ f(x, y) = (1/2) \iniy^ x [f(x, y) + 
f(y, x)] is the symmetric limit. Clearly we can multiply any fermion bilinear by any func- 
tion F(x,y) such that F(x,x) = 1. As usual, F(x,y) can be determined by demanding 
that the different observables obey physical conditions such as current conservation or 
gauge invariance [ 34| . Consider for example the current j^x), which is the relevant one 
here, and suppose we just choose F(x,y) = 1. Recall that the relevant theory here, ac- 
cording to ( p,18| ) for p = 0, is that of a fermion field in the presence of an external vector 
field A^x). Then classically, j M (x) is a conserved current in such a theory. However, it is 
easy to check that <9£lim* ^(y)^(x) = ]im y ^ x (l/2)A fM (x)[^(y)j fl 'i(;(x)-^(x)j tl 'ij;(y)], 



i.e, that choice for F(x, y) would not lead to a classically conserved current, 
instead that the correct regularised current is 



lim s ?/%)7 M ^(a;) exp 



ig I A a (Qd£ a 
y 



We see 



(4.10) 
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which is classically conserved before taking the symmetric limit. Also ( |4.10| ) can be seen 
to satisfy the Ward identity <C j T u e9 (x) 3>= 0, just by writing -C ^ a (x)^p(y) ^>= 
G>i3(x,y), the advanced exact fermion propagator, and, in addition, when dealing with 
a gauge theory, this definition of the regularised current is gauge invariant, another re- 
quirement the current should satisfy |L9|, |25], j34[ . 



We note that in the free case g = 0, the point-split current yields automatically a 
finite result for correlators with current insertions, without any need of renormalisation. 
The reason is that, from ( |3.6| ), we see that the divergent part of S± always behaves like 
(e ± ex) -1 with e u = y u — x^ and hence it cancels when the symmetric limit is taken. 
Furthermore, we will now see that those are the only infinities also for g ^ so that no 
extra renormalisations are needed and the choice ( |4.10| ) for F(x, y) is the only one giving 
the right equivalence with the SG model. 



4.3.2 Chiral invariant mixed correlators 

The Thirring model Lagrangian and the fermion measure are invariant under charge 
conjugation C and therefore <C ^>mt = 0, consistently with e nu d u <C <p(x) ^sg= 0, 
following just from the <p — > —0 invariance of the SG model. Thus, the equivalence is 
too simple with just one insertion of the current. In this section we concentrate on the 
mixed propagator with one current and one cr + cr_ insertion. This correlator again has 
the useful property that it is chiral invariant so that the selection rules are very simple. 
Notice that our results can then be extended without difficulty to an arbitrary string of 
(j + o"_ and/or j^'s or even to chiral non-invariant combinations, but for clarity we restrict 
here to the simplest example. In turn, we will be showing that the equivalence is also 
true for mixed mass-current correlators. 

Thus, according to the discussion in the previous section, consider the correlator 

F ix (x,zi,z 2 ) =<C T c f^ 9 (x)(j + (z l )a_(z2) >mt 



with j r a e9 (x) in ( |4.1C ). After similar manipulations as those of the previous sections one 
obtains 



F fl (x,z 1 , z 2 ) 



Z F (T) 



00 ml n 
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n+1 



5^(x 1 -z 1 )6^(y l -z 2 ) 
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rx o n + 1 

< T c e TK[<t,{x) ~' , ' {y)] e iKea/3 -'y 9 -q e -M<K*i)-4><3ii)] > c 

3=1 



(4.11) 



where the rj field dependence is the same as in ( |3.21| ) and has been integrated out. Now 
notice that the only pole in (x — y)^ = e a — > is in the fermion correlator in ( |4.11 ) 
and that it goes like e _1 . Thus the first two exponentials in the boson free correlator in 
the last line can be expanded up to 0(e). Let us analyse first the fermion free thermal 
correlator in (|4.11| ). After shifting ip — > ipj , all the required fermion correlators are of 
the form 



(_l)n+l Tc 



n+1 



3=1 



26 



23 



2n+3 
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y.i 
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0{e 2 ) (4.12) 



Q 



where Zq are defined in (|2.19|) , we have used ([3 .71) and (|3.8|) and we have expanded to 
the relevant order in e. As for the boson correlators, to 0(e), we only need to calculate 



n+l 



< Tcdfflx) J] exp [-2k(0(^-) - </>(%•))] > = z z; ( /r/f () - ; 
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where 

J K (z) = -2i £ [*<*>(* - x,) - 5^(z - %•) . 
i=i 

Notice that by taking the partial derivative we make sure that there is no extra jl- 
dependence coming from the propagator. Also recall that the UV e factor above (not 
to be confused with e M ) will be absorbed in the renormalisation of the mass and the a± 
fields as indicated in previous sections. Finally, use of ( [2.18|) in (|4.12|) gives 
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(4.13) 



with < T c cxf (21)0-^(^2) > m = o in (O) for N = 1. Notice that, thanks to ( gig ), the 

above result is independent of the way we take the limits eo — > and ei — > 0, which is a 
good check of consistency. 

On the other hand, in the SG model we have to evaluate 
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Q 2 (x - 2 2 ) ^ Q 2 (x - j/,) 

with J\{z) = —i\(5^(z — 2^) — 5^ 2 \z — z 2 )) and Af(z±, z 2 ) as in (|4.7|) . Thus comparison 
of ( |4.14|) with (|4.13j) and use of the equivalences we have already proved gives 

« T c j r ^(x)a R (z 1 )a R (z 2 ) » MT = {£) ^ « T c d^{x)A R { Zl )A R {^) »sg 
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which proves the current equivalence for this correlator. 



5 Bosonisation of the massive Thirring model at 
nonzero chemical potential 

In this section we study the thermal bosonisation of the Thirring model at nonzero 
chemical potential fi. The chemical potential is the Lagrange multiplier associated with 
the conservation of the number of fermions minus anti-fermions (net fermion number) 
Qf = J dxj°(t, x). Hence, the Lagrangian density of the /i ^ Thirring model is just 



(5.1) 



By virtue of the current equivalence we have proved in section [4.3| , we can say that at 
/i = 0, < Q F » >a/t =< Qk* >5G, where 



Q 



K 



^^^^ 



2^ 



[<f){t, +oo) - <j){t, -oo)] 



corresponds classically to the conserved number of kinks minus anti-kinks ||. Recall that 
classically <fid(t,±oa) = 2nn±/\, with n± G 2, so that Qk[4>ci] = n + — n -- We regard 
this equivalence as the thermal identification between Thirring fermions and SG kinks at 
T / 0, thus extending the T = standard result |], [|, discussed in the introduction. 
Physically, we interpret it as the existence of fermion excitations (the kinks) in the thermal 
bath. 

The presence of excitations with net fermion number becomes crucial when chemical 
potentials are included and hence the grand-canonical ensemble partition function is 
considered. As the averaged net fermion number density p(/i) = (/5L) _1 C Qf > will 
in general be different from zero, a natural question to ask within this context is what 
is the bosonised version of the theory for /i / 0. As we commented in the introduction, 
the answer has been given in |19| for the massless case, where the free boson partition 
function acquires an extra //-dependent term, unlike, for instance, the massless Schwinger 
model, where there are no elementary fermion excitations in the bosonised theory (all 
the electric charge is confined), even for nonzero T and fi |19]. The Schwinger model 



partition function is then equivalent to that of a boson field with mass e/y/ir, with e the 
electric charge, and the chemical potential dependence vanishes in the partition function, 
yielding p(/i) = 0. 

Our objective in this section is then to analyse the partition function of the \i ^ 
massive Thirring model ( |5.1| ), 



dipdip exp 
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d 2 x£ M T[ip,ip] fA 



(5.2) 



The simplest way to proceed is to make use of the current equivalence proved in section 
|4~3| . Formally the argument goes as follows: we expand Zmt(T,(i) in /x, so that every 
term in the expansion is a \x = correlator of a product of j°{xj) at different space-time 
points Xj. Then assume that the T^O current equivalence ( |1.4|) holds for any arbitrary 
product of currents so that we conclude 
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with 

Zsg^(T,h) = Np [ d<fi exp - f d 2 x ( C SG ~ Tr^-T^x) 

J periodic JT \ 2lX OX 



(5.4) 



/periodic 

i.e, the bosonised action is the SG model plus an extra term which is topological, in 
the sense that it only depends on the value of the field at the spatial boundary (x = 
±00). As we have commented before, that term is interpreted as the result of excitations 
with net kink (fermion) charge being present in the thermal bath and having associated 
the chemical potential \x in the grand-canonical ensemble. Recall that an analogous 
contribution was found in the chiral Lagrangian for low-energy QCD [IH] where the 
role of kinks (Thirring fermions) is played by the skyrmions (QCD baryons): the chiral 
Lagrangian at finite baryon density acquires, among other things, a new factor fiQsx 
with Qsk the skyrmion topological charge, analogously to ( |5.4|) . The existence of this 
term provides an alternative proof, based on Statistical Mechanics, of the old result that 
the classical kinks have their quantum counterpart in a particle with fermion number 
equal to one jl], and it is indeed another example of /^-dependent topological terms 
found in the literature fT7|j . 

The equivalence (|5.3| )- (|5.4| ) is the central result of this section. We want to stress that, 
formally, we have proved it only in perturbation theory around the m = a = theory. 
Indeed, we can perform a simple consistency check of ( |5.3|) , since for a = 0, Z SGlL should 
give the answer for the massless Thirring model obtained in ]HJ. In that case, JCgQ in 
( |5.4|) is just a free boson Lagrangian, so let 

4>{x,t) ^<f>(x,t) + ^x. (5.5) 

As for the field measure, in the free theory the field is allowed to have A<ft = <p(L/2, t) — 
0(— L/2,t) = c with c 6 R arbitrary and L the length of the system, so that the total 
energy still remains finite when L — > 00. In other words, the free theory has trivial 
topology since any vacuum configuration can be reached by infinitesimal variations of 
the field. Therefore, we can write the path integral measure formally as J R dc JacA=c d<j)- 
Thus, under (|5\5|), A0 — > A(f> + fi\L/2ir, but then, by shifting c — > c — fi\L/2n, the path 
integral measure remains invariant. Therefore, by completing the square and using ( [L.5|) , 
we find 



(3L- ^ 



2(<K + g> 



Z B (T) (5.6) 



which is exactly the result found in |19| for the massless Thirring model partition function. 

In what follows, we provide an alternative way of obtaining the result ( |5.4j ) without 
having to expand in \x. For that purpose, we will make use of some of the results obtained 



in [y| for the massless case. Hence, let us return to the massive theory for /i ^ and 
compare the series expansion of (|5.2j ) in m with that of ( |5.4j ) in a, in the same way as 
we did in previous sections for \i = 0. Notice that we expand now around the massless 
fi 7^ theory, without expanding in /x. An important remark is in order here: for every 
term in the a-expansion, which is an a = boson correlator, we will consider that 
the measure is translation invariant, as discussed before. In the mass expansion of the 
/i ^ Thirring model, we will make the same assumption whenever a boson correlator 
appears. Therefore, it should become clear that we are perturbing around both the a = 
Lagrangian and path-integral measure and hence our results are valid only in perturbation 
theory. It remains unclear to us whether we would arrive to the same conclusion, namely 
(|5.4|), within a semi-classical approximation, i.e, perturbing now around a nontrivial 
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classical configuration, which would be consistent for A <C 1 |], [22| but that is beyond 
the scope of this work and it will be analysed elsewhere. 

Therefore, as in previous sections, let us introduce the auxiliary field A a and expand 
in m in ( |5.2|) . Then we have 

Z j\(T^ Jperiodii 
oo ^,2n 



Z M T{T,fi) 



f dA^exp -- [ d 2 xA a (x)A a (x) Z f [A;T,fi} (5.7) 

J periodic \_ Z JT i 

Z f [A;T,fi] = J2r~J\2 / d 2 x j d 2 y j a + (x j )(T-(y j ) 

n [fT'-l Jantiper j_-^JT 

d 2 xip ($ + ig 4(x) - /ry°) 



n=0 

x exp 



(5.8) 



Notice that (|5.7| ) and ( |5.8|) are nothing but the mass expansions of ( [3.17| ) and ( |3.18|) 
respectively for fj — r] — 0. In order to calculate the m = correlator in Q5.8| ), we put 
back the external sources f] and rj and take mo = in (|3.18|) and recall the following 
result in JljJ 



Z7 = [Ar/,r/;T,/i] 



exp 



r 



cPx / d 2 yri(x)G(x,y,A;[i)7](y) 



T 



x exp 
x exp 



~2 //^ j T *vM*)K*l? - v)Mv) 

—i— d 2 xA Q (x 

IX JT 



(5.9) 



Here Zq(T,[i) is the /z 7^ free massless fermion partition function which in 2D is given 
by 



logZ F (T,/i)=logZ F (T,0)+/3L 



2tt 



and 



PaP/3 
p 2 



with p Q = 2irnT, 116Z. G(x, y, A, /i) is the exact massless fermion propagator, which is 
given by \T§\ 

G(x,y,A;fj,) = exp {ig[x(x) - x(y)]} S(x, y; fi) 



T 



d 2 zA T (x - z) $ z 4(z) 



with A? the free boson propagator in ( p,14| ) and S(x,y;fi) is the free fermion /1 / 
propagator, solution of (^ — yU7°)S'(a;, y; /1) = 5^ 2 '(x — y), for which we can write 



S(x,y;n) = exp \-i{x -y J S(x,y;0). 

Now redefine fields and sources as ip ~ > "07° an d f] — > ^7°, so that 770 — > 771 ^2 +^2^1 an d 
■077 — >■ 0i772 + 02771 and, with the above definitions, 



r](x)G(x,y,A;fi)r](y) = J2Va(x)G a (x,y } A; n)rj a (y), 

1 

1 ZZ(x-y) 



a=l 



G±(x,y,A;fi) 
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exp 



X ± (x)-X ± (y)Tt^x L -y i ),(5A0) 



where G + = G\, G_ = G2, Zq are given in fl2.19|) and 



ig [ d 2 zA T (x - z) \H(z) ± iE(z)\ 

JT 



with H = doA + d\Ai and E = doAi — diA . Now, let us decompose the vector field as 
in ( |3.20| ), but letting now A0(x) 7^ 0, (see our previous discussion for the massless case), 
so that we have 

~\L d2x St rf Vu^it^x - y)Mv) = ^ I d2xd *<pd a <P 

d 2 x [d a r]d a 7] - d a (j)d a <j)} 



-- I d 2 xA a {x)A a {x) 



2Jt 



—1- 



.1*9 



7T JT 



d xAn(x) 



2g 2 Jt 

/i 



d 



g J d 2 zA T (x - z)E(z) 



d2x a i 

7T JT OX 

-id>(x) 



[X 



(5.11) 



and H will not play any role. 

The final step in order to calculate the cr + <7_ correlator in (|5.8|) , is to differentiate 
( [5.9D with respect to the sources fj a and rj a , after redefining both fields and sources as 
indicated above, that is, 



<<T c n^>-(%)>A 
3=1 



-> « T c ^j)M^)My^2(yj) »7 =0 



where the subscript A indicates the thermal expectation value in the presence of A^, 
according to Q5.8|) . Then, taking into account (|5.10|) and (|5.11j) , we find 
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periodic 



exp 



3=1 
x exp 



n „ 

JJ y d 2 Xjd 2 yj 



d 2 x yl- 2 d a <P{x)d a <P{x) - ~<f>{x) 
Uk<j [Q 2 (xj - x k )Q 2 (yj - y k )} 

nu [g 2 (x,-i/ fc )] 



{2 ? [0(^.)-0( % .)- / ,( a; )-^)]} 



(5.13) 



At this point, and to solve the path integral in <p above, we make use of the measure 
translation invariance argument discussed before and, accordingly, we shift — > in{(p — 
ifinx 1 /it) so that the /1 term is eliminated in the boson action in (|5.13|) . Then, we can 
calculate the boson correlator as in previous sections and finally obtain 



Z MT (T )M ) = Z^°(T^)Y: 



n=0 



ni 



m / p 
2/3 \ f 



K 2 /n 



2n 



3=1 t 



X 



exp 2i/i ^1 - ^ (xj - yj) 



(5.14) 



31 



where Z^}jP{T,fi) is given by (|5.6|) (recall that Zq(T) = Zq{T)). We readily check that 
the above yields ( |3.22| ) for fi = 0. Another consistency check is that the fi 7^ partition 
function we obtain is real, which can be easily checked in ( |5.14| ) just by redefining inte- 
gration variables Xj <-> yj. Besides, by taking only the lowest terms in the m-expansion, 
we have numerically checked that it is also positive definite. 

Now, we expand (|5.4|) in «o and perform the shift ( |5.5| ), assuming the measure invari- 
ance. Then, using ( p.6|) and (|1.5|) - (|1.6|) , we obtain exactly ( |5.14|) and then the equivalence 
rSD-fFl follows. 



6 Conclusions and Outlook 

We have studied in detail bosonisation at finite temperature and nonzero fermion chemical 
potential in the massive Thirring and sine-Gordon models. The main new results of this 
paper are the following: the extension of canonical methods to T > and fi = to 
prove that the partition functions of the two models are the same; the proof, via path 
integral methods of the equivalence between thermal averages of different sets of zero- 
charge operators in the two models; and the analysis of the fi 7^ case where the main 
result is a topological new term in the bosonised SG partition function, accounting for 
the existence of fermion modes in the spectrum. All our results have been proved in 
perturbation theory around the massless Thirring model or, equivalently, the free scalar 
theory. 

In the canonical operator proof, the use of the thermal normal ordering was crucial. 
Also our results for the operators (mass and current) equivalences, as well as for the 
T > 0, 7^ partition function, clearly show the existence of fermion excitations in the 
thermal bath. This constitutes the thermal identification between sine-Gordon kinks and 
Thirring fermions, so generalising the existing knowledge about this subject obtained 
when Statistical Mechanics features were not considered. In addition, the topological 
term obtained for fi 7^ is the natural analogue of a similar term found in the effective 
low-energy action for QCD at finite baryon density. 

Among the issues which we would like to investigate further in the future is the 
semiclassical approach, which we have commented about in the text, especially in what 
concerns the /i 7^ case and the topological action. Also, as was noted in the introduction, 
we would like to exploit the relationship between fermions and kinks, when these are 
treated as topological defects, since we believe that these 2D systems could become a 
very promising playground for testing cosmological models. The experimental realisation 
on Josephson junctions could then provide precise laboratory tests where, in particular, 
the results obtained in this paper could be physically probed. Other points of interest 
which which we are looking at in more detail at the moment include a numerical analysis 
of the free energy and the fi 7^ fermion density, the high temperature limit of these 
systems, and their behaviour near the point A 2 = 8tt. 
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